home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / clahqr.z / clahqr
Text File  |  1996-03-14  |  4KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCLLLLAAAAHHHHQQQQRRRR((((3333FFFF))))                                                          CCCCLLLLAAAAHHHHQQQQRRRR((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CLAHQR - i an auxiliary routine called by CHSEQR to update the
  10.      eigenvalues and Schur decomposition already computed by CHSEQR, by
  11.      dealing with the Hessenberg submatrix in rows and columns ILO to IHI
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE CLAHQR( WANTT, WANTZ, N, ILO, IHI, H, LDH, W, ILOZ, IHIZ, Z,
  15.                         LDZ, INFO )
  16.  
  17.          LOGICAL        WANTT, WANTZ
  18.  
  19.          INTEGER        IHI, IHIZ, ILO, ILOZ, INFO, LDH, LDZ, N
  20.  
  21.          COMPLEX        H( LDH, * ), W( * ), Z( LDZ, * )
  22.  
  23. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  24.      CLAHQR is an auxiliary routine called by CHSEQR to update the eigenvalues
  25.      and Schur decomposition already computed by CHSEQR, by dealing with the
  26.      Hessenberg submatrix in rows and columns ILO to IHI.
  27.  
  28.  
  29. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  30.      WANTT   (input) LOGICAL
  31.              = .TRUE. : the full Schur form T is required;
  32.              = .FALSE.: only eigenvalues are required.
  33.  
  34.      WANTZ   (input) LOGICAL
  35.              = .TRUE. : the matrix of Schur vectors Z is required;
  36.              = .FALSE.: Schur vectors are not required.
  37.  
  38.      N       (input) INTEGER
  39.              The order of the matrix H.  N >= 0.
  40.  
  41.      ILO     (input) INTEGER
  42.              IHI     (input) INTEGER It is assumed that H is already upper
  43.              triangular in rows and columns IHI+1:N, and that H(ILO,ILO-1) = 0
  44.              (unless ILO = 1).  CLAHQR works primarily with the Hessenberg
  45.              submatrix in rows and columns ILO to IHI, but applies
  46.              transformations to all of H if WANTT is .TRUE..  1 <= ILO <=
  47.              max(1,IHI); IHI <= N.
  48.  
  49.      H       (input/output) COMPLEX array, dimension (LDH,N)
  50.              On entry, the upper Hessenberg matrix H.  On exit, if WANTT is
  51.              .TRUE., H is upper triangular in rows and columns ILO:IHI, with
  52.              any 2-by-2 diagonal blocks in standard form. If WANTT is .FALSE.,
  53.              the contents of H are unspecified on exit.
  54.  
  55.      LDH     (input) INTEGER
  56.              The leading dimension of the array H. LDH >= max(1,N).
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCLLLLAAAAHHHHQQQQRRRR((((3333FFFF))))                                                          CCCCLLLLAAAAHHHHQQQQRRRR((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      W       (output) COMPLEX array, dimension (N)
  75.              The computed eigenvalues ILO to IHI are stored in the
  76.              corresponding elements of W. If WANTT is .TRUE., the eigenvalues
  77.              are stored in the same order as on the diagonal of the Schur form
  78.              returned in H, with W(i) = H(i,i).
  79.  
  80.      ILOZ    (input) INTEGER
  81.              IHIZ    (input) INTEGER Specify the rows of Z to which
  82.              transformations must be applied if WANTZ is .TRUE..  1 <= ILOZ <=
  83.              ILO; IHI <= IHIZ <= N.
  84.  
  85.      Z       (input/output) COMPLEX array, dimension (LDZ,N)
  86.              If WANTZ is .TRUE., on entry Z must contain the current matrix Z
  87.              of transformations accumulated by CHSEQR, and on exit Z has been
  88.              updated; transformations are applied only to the submatrix
  89.              Z(ILOZ:IHIZ,ILO:IHI).  If WANTZ is .FALSE., Z is not referenced.
  90.  
  91.      LDZ     (input) INTEGER
  92.              The leading dimension of the array Z. LDZ >= max(1,N).
  93.  
  94.      INFO    (output) INTEGER
  95.              = 0: successful exit
  96.              > 0: if INFO = i, CLAHQR failed to compute all the eigenvalues
  97.              ILO to IHI in a total of 30*(IHI-ILO+1) iterations; elements
  98.              i+1:ihi of W contain those eigenvalues which have been
  99.              successfully computed.
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.